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Abstract 

Criteria are given by which dissipative evolution can transfer populations and coherences between 
quantum subspaces, without a loss of coherence. This results in a form of quantum error correction that 
is implemented by the joint evolution of a system and a cold bath. It requires no external intervention 
and, in principal, no ancilla. An example of a system that protects a qubit against spin-flip errors is 
proposed. It consists of three spin 1/2 magnetic particles, and three modes of a resonator. The qubit 
is the triple quantum coherence of the spins, and the photons act as ancilla. This article is a greatly 
expanded version of a letter submitted to Physical Review Letters. 

PACS number(s): 03.67.Lx 

1 Introduction 

Quantum computation is of interest because algorithms have been discovered with a significant speed-up 
over any classical algorithm J3J, although these may be unique cases ||. It is very likely that any 
physical implementation of a quantum computation will require some form of active quantum error correction. 
Quantum error correcting codes (QECC) have been devised [0, |[ @ and experimentally demonstrated 
H that can protect a set of states, the codewords, against a set of errors. QECC is similar in spirit to 
quantum erasure experiments ||, but with the twist that one is not allowed to manipulate the environment. 
The surprising fact is that one can still disentangle the codewords from the environment, by transferring the 
entanglement to another set of states, the ancilla. 

However, implementing QECC is a formidable task. There is a high premium placed on using as few 
qubits (two-level systems) as possible, because as quantum systems grow in size, the number of transitions 
to be manipulated, unwanted thermal effects |fTo|| , and decoherence rates jllj all increase exponentially. But, 
to take a specific example, the fault-tolerant error detection and repair of even a single qubit can require 15 
physical qubits, 5 to store the two codewords, and 10 of which must be in known states of zero entropy ||. 
In addition, depending upon how one counts a "logic gate" , as many as 28 coherent manipulations of pairs 
of qubits are required for each repair, because "measuring the stabilizer" means finding the eigenvalues of 
operators such as I x ilx2lz3lz5 (see Fig. (2) of Ref. [|)). Such control over a 32,768-level system is a daunting 
task, even for a highly coherent spectroscopy such as NMR. Although the efficiency of QECC improves for 
larger computations, a physical scale- up factor of 22 is still required to factorize a thousand-digit number 
[ p2[ . Part of the difficulty stems from the need to know which error has struck, in order to repair it. This 
is because different errors rotate the codeword states about separate axes in Hilbert space. By containing 
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information about which error occurred, the ancilla also provide a conditional axis about which rotation can 
coherently repair an error. Although this seems like an air-tight argument, there is another way to approach 
QECC, which we will explore here. 

To begin with, note how curious it is that QECC can assign a unique status to the codewords. While a 
classical probability space inherently contains a privileged basis, Hilbert space does not, and this difference 
has some striking consequences jf3|. In order to function, QECC requires access to ancilla in a state of 
zero entropy K3, which suggests that one could view QECC as a controlled cooling of the system. It is 
dissipative evolution that adds classical aspects back into Hilbert space. A large body of work exists that 
model a diverse range of relaxation phenomena in magnetic Jl5| and optical resonance [fl6[ 17| . They use 
Lindblad equations of motion |ll[| . An earlier approach that did use a Lindblad equation fl8 ] implemented 
QECC as a limit of very fast external manipulation. In contrast, we seek an approach that is distinct from 
the concepts of error detection and repair. 

We show here that dissipation can be used to implement an "automatic quantum error correction" 
(AQEC), so called because error correction results exclusively from the joint evolution of a system coupled 
to a cold, Markovian bath. No intervention is required by the programmer, and in theory, no ancilla are 
required, although this would be unlikely in practice. Clearly, dissipation can be used to stabilize two distinct 
states of a quantum system that could store a classical bit of information. What is not obvious, is whether 
such a system could also hold a qubit, since dissipation usually destroys coherence. The key ideas are to use 
codewords such that errors must add energy to the system, and to set up the evolution of the system such 
that excitation and environmental entanglements are expelled from distinct codewords in a symmetric way. 
This prevents the bath from gaining information on the codewords, and thus coherence can be maintained. 
In the last section, we outline a system that utilizes three magnetic spin 1/2 particles, and three photons, to 
implement an AQEC that protects against spin- flip errors. It requires only well understood interactions from 
magnetic resonance spectroscopy, and is intended to show that AQEC has potentially real- world applications. 



2 A simple QECC example. 

To begin, let us review the idea of quantum error correction by way of a simple method that protects a 
single quantum state against environmental entanglements [|lQf| . The idea is similar to that of a quantum 
eraser experiment J^], but with the twist that one is not allowed to interact the environment. Three two- level 
systems, labeled as S, A and E, are initially in the state (a\ls) + &|0s))|0£:) \0a)- The goal is to keep S in 
its current state. An interaction between S and E creates the new state (ci\1s}\pe} + 6|0s) |?b))|0a)- The 
environment is scattered into two states, \pe) and When \(pE\qE)\ < 1, the final state of E depends 

upon the initial state of 5 1 , so they are entangled. If \pe) = — |<Zs), the phase of S has been flipped. To 
repair S, first note that the entangled state can be written as: 

&|0 S >) (\pe) + |sb>) + falls) - &|0s>) (\pe) - 

Suppose we can externally manipulate the qubits. Conditionally flip A, if the sign of the state S is flipped 
from what we expect it to be. In the language of QECC, this is "measuring the stabilizer", or "detecting 
the error" . A serves as the memory. Next, flip the sign of S, conditional on if A detected an error. This 
is "repairing the state" . Both of these actions are unitary transforms on S and A only; the environment is 
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never directly manipulated. After measuring the error, the new state is: 




b\0s) \pe) 




and then, after the repair, 



alls) + &|0s>) \\ (\pe) + |te>) |0. 4 > + (\pe) - |«b)) |U> J 

The original state of 5 has re-emerged! The entanglement between S and i5 was transferred to be between 
A and E, without ever touching E. In order for this scheme to work, it is crucial that A is initially in a 
single pure state, or in a state of zero entropy. We expect that we can achieve this by cooling A down to 
°K by the third law of thermodynamics. However, there are systems such as protons in ice or frustrated 
spin lattices |2(J that are postulated to violate the third law. Since cooling these systems still leaves them 
in a state of non-zero entropy, one should avoid using them as ancilla. 



2.1 A dynamical re-formulation of the above example. 

The next step is to transform the above error correcting method into the language of a quantum system, 
relaxing towards equilibrium. We will make use of the operator formalism of NMR pl| . The qubit states 
are |0) and |1), and the projection operators are I a = |1)(1| and Ip — |0)(0|. The raising and lowering 
operators are I + = |1)(0| and 7_ = |0)(1|, respectively. The Hermitian Pauli operators are I x = (i+ +J_)/2, 
I y = (J+ — i_)/2i, and I z = (I a — Ip)/2, and / is the vector formed by them. The subscript also indicates 
which spin is acted upon, so I n x acts only on spin n. This section is similar to that of Ref. |l8| ], but with 
the difference that the measurement of the syndrome and the repair process are treated more explicitly. 

The Hamiltonian of Eq. (|l|). is designed to continuously implement the example of the last section. To 
keep S in the state first flip A at the rate d, if S departs from \ls)- If A has flipped, then S is flipped at 
a rate r. To complete the process, A is cooled at a rate c by interaction with a bath of harmonic oscillators 
with a broad spectral response. If the bath temperature is low in comparison to the separation of the levels 
of A, then the density matrix p evolves as JO], psf , 

H = r{lA,0 + iA.als.x) + d(I S .a + Is.plA,x) , 

d ( ^ ) 

= -i[H,p] - c(I AtCt p + pI A>a - 2I A _ pI A ,+) 

We suppose the errors occur rapidly in comparison to the system dyanmics, so they are modeled as instan- 
taneous transforms. However, slower interactions can also be corrected Q. 

Dissipative evolution is usually handled in Liouville space, where p is a vector, and transformations like 
—i[H,p] are matrix- vector multiplications p2[ . These matrices are called superoperators, since they operate 
on operators. Eq. (|l|) becomes a set of linear differential equations, p — Tp. The elements of T are then 
indexed by how they transform the populations and coherences of an orthonormal set that spans the Hilbert 
space: each element of T transforms a \j)(k\ to a |n)(m|. When c = 0, the evolution is unitary, and T has 
two kinds of eigenvalues: A = 0, corresponding to populations of eigenstates of H, \n)(n\, and A = ±i/3, 
corresponding to coherences |n)(m| and |m)(n|. Unfortunately, Liouville space also increases the problem 
size: N qubits now require an evolution superoperator with A N eigenstates. 

When evolution is dissipative, T is not a symmetric matrix. It can still be written as the outer product 
of its right and left eigenvectors, T — Y] X n f n <g> l n , but in general the r n are not orthogonal. However, 
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In ' — 5nm, which allows one to formally solve the equation of motion for an operator as x(t) = r n (l n ■ 
x(0) ) exp(A„t). The structure of Eq. ([!]) implies that T conserves tr(p(t)), but a pure state will not 
necessarily remain pure. 

Does the system of Eq. (Q) work? Fig. ([!]) plots the real parts of the eigenvalues of T of Eq. (|l|) for 
various values of d, r and c. The only stable states of T will be those with 3?e(A n ) = 0, and there is only one 
such state: |1s0a)(1s0^|. This is how dissipative evolution can confer a privileged status on a state. 

A numerical integration of the system dynamics also shows this. Fig. (Q) plots the linear entropy, defined 
as < tr(p(t) — p 2 {t)) < 1, of p(t), starting from the corrupted state |0s0a) (0s0a|- The rate at which 
the error is repaired is dominated by the eigenvalue of T with the least negative, but nonzero, real part. 
Curiously, a larger c is counterproductive, as it traps the state into a cycle: 

\0 s 0a) d6to l7 1°^) repairs \1 S 1 A ) cool -» |1 S A ) 
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Fig. |l|. The real parts of the 16 eigenvalues of 
the superoperator T for some values of d, r and c. 
Note that there is only a single stable state. 
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Fig. |2|. The linear entropy during the continuous 
error correction, for d, r, c = 1, and doubling each 
parameter separately. The starting state is p = 
|0 s 0a)(0 s 0a|. 



3 Conditions for AQEC. 
3.1 Repairing the Populations. 

We now show what conditions are necessary in order that dissipative evolution can automatically protect 
a subspace of codewords against a given set of errors. We first suppose that the system obeys a Lindblad 
equation of motion. In general, this is not a trivial assumption pT[ . The most speculative condition in 
deriving a Lindblad equation is that the system and the bath initially factorize. Curiously, it can be justified 
here on the grounds that a properly working error correction should drive the system to this state. A more 
troublesome condition is that if degenerate transitions are coupled to the Markov bath, they must couple to 
orthogonal bath modes . 

Two conditions can be stated immediately. We arc assuming that evolution for a sufficient time, T, can 
repair any error. Thus, exp(r T) must be a repair superoperator. Necessary and sufficient conditions for 
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its existence are known from QECC |iS| ]. Under these conditions, the original codeword populations and 
coherences might be transported elsewhere in Hilbert space, but they are not destroyed. The second condition 
is that the codewords must be immune from the influence of the bath, or that they form a decoherence free 
subspace with respect to the system / bath coupling |^5| . 

An example is instructive. Suppose that we are interested in protecting a two-codeword system against 
spin-flip errors. The system is split into two groups of qubits, S and A, where the A are continuously cooled. 
This is not necessary for the general argument, which can be formulated entirely in terms of the eigenvalues 
of r. However, it simplifies the physical interpretation. The system evolution is given by |l6|, [l7j 

dp ancilla , v 

-^=-i[H,p]- ^ Cn[ln,aP + pIn,a-2In,-pIn,+ ) ( 2 ) 

where H acts on both the S and A. The second term irreversibly draws population from the |1) states of 
the ancilla spins, and places it in the |0) states. Choosing codewords of the form \iI>ii)\0a), for which (1) 
the A are in their ground states, and (2) the codewords are eigenstates of H, will satisfy the criteria for the 
decoherence-free subspace. 

The need for the QECC conditions can be seen as follows. Suppose we choose the two-S* states |00) and 
|11) as the codewords. But then the errors Ii jX \0Q) and 22,ic|ll) both result in the same state, |10). Under 
the Markov approximation, the system can not know which codeword was the original codeword, and so T 
can not repair these errors. But the three-S states 1 000) and 1 1 11 ) will work, since the spaces spanned by all 
the errors acting on each codeword are now disjoint. Thus, errors should transfer separate codewords into 
disjoint subspaces. 

Now consider how to repair the codeword populations. The set of errors acting on a codeword, and all the 
further states that the corrupted codeword evolves into under T, form a subspace, as indicated in Fig. (||). 
Call this subspace the "funnel" associated with the codeword, but excluding the codeword state itself. The 
name is suggestive of its role in AQEC. The QECC conditions already require the initially excited states to 
be disjoint between separate codewords. Thus, if we add the third condition that T draws all population 
from each funnel state into its associated codeword, and transfers no amplitude between funnels, then the 
codeword populations are repaired. 




Fig. |3|. The level diagram of a hypothetical sys- 
tem. The errors (solid arrows) transfer amplitude 
into the disjoint funnels associated with each code- 
word. The three levels between the brackets form 
the funnel for the labeled codeword to the right. 
Dissipative cooling of selected transitions (dashed 
lines) then returns the populations to their origi- 
nal codewords. 



codeword 



An important difference between the usual method by which QECC is implemented, and AQEC, has 
emerged. By placing the burden of the repair on the system / bath coupling, AQEC, in theory, requires 
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no ancilla. Consider how to repair the error Ii^ x , acting on the three- S and two- A codewords 1 000, 00) and 
| 111, 00). Suppose that the unnormalized states |100, 00) ± |000, 10) and |011, 00) ± | 111, 10) are eigenstates 
of H. The corrupted state 1 100, 00) now periodically becomes the state |000, 10), where cooling of the 
first ancilla returns it to the codeword |000, 00) . However, we can not choose |010, 00) ± |000, 10) and 
1 101, 00) ± 1 1 1 1 , 10) as eigenstates of H in order to repair the error I 2 , x , because they are not orthogonal to 
the first set. It seems that we must choose |010, 00) ± |000, 01) and 1 101, 00) ± | 111, 01) to repair I 2 , x , and 
|001, 00) ± |000, 11) and |110, 00) ± 1 1 1 1 , 11) to repair I 3yX . In this strategy, we must be able to distinguish 
between the different errors in order to be able to repair them, and counting e different errors requires log2£ 
binary digits (or ancilla qubits). 

But the third condition for AQEC only requires that H mix the states |100, 00), |010, 00), and |001, 00) 
with |000, 10), even if only partially! Because cooling irreversibly draws probably away from the excited 
ancilla, any degree of mixing will do. In other words, when the 4x4 block of H corresponding to the above 
states is diagonalized, each eigenstate should have a non-zero projection onto the state |000, 10), and similarly 
for the other codeword. We still require at least one ancilla here, because the system was split into S and A, 
and only the A are cooled. If a bath / system coupling is found that directly cools the funnel to codeword 
transitions as in Fig. (||), then no ancilla are necessary. However, not all the conditions for AQEC have been 
stated yet. The rest of these come from the seemingly bizarre notion that we can use dissipation to restore 
a coherence. 



3.2 Repairing the Coherences. 

The QECC conditions ensure that the codeword coherences are transferred, but not "measured" , by the 
environment. AQEC must transfer them back. What happens to coherences during dissipative evolution is 
a subtle point, which is best explored by way of a comprehensive example. Codewords of the form \ip n ) 1 00) , 
with two ancilla, will serve this purpose. The environment, |e), is initially unentangled with the computer. 
An interaction, U, can entangle the system so that p6[ 
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After the error, the amplitude originally in each codeword is spread throughout its funnel. While the funnel 
states {|</>n^)} can be chosen as an orthogonal set for each n, this is not true in general for the {|e^)}. 

As yet, there is no constraint on either how separate codewords can excite the ancilla qubits, or how 
the dynamics of the repair should proceed. Suppose H uses the first ancilla to repair the codewords n = 
and 1, the second ancilla to repair n=2, and both ancilla to repair n=3. That is, H mixes each |4 fe) )|00) 
with |^o)|10), and so on. Let us follow an argument analogous to the "quantum jump" approach |[L7| . 
The relaxation process is divided up into small time steps, At, during which the system and bath evolve 
separately. At the end of each interval, a fraction of the amplitude in each excited ancilla state jumps into 
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a de-excited state. Tagging on two more qubits to represent two modes of the cold bath, at the end of each 
time interval, a fraction of the amplitudes make the following jumps: 
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The entire process is repeated until a time T, when the relaxation process is complete. 

The heart of the argument relies on the idea that, in the limit of a large number of cold bath modes 
interacting with the ancilla, it is very likely that different ancilla that de-excite at different times, will transfer 
their excitation to orthogonal modes of the bath. Once excited, these modes do not further influence the 
evolution of the computer, i.e. there is no back- reaction from the bath. In this case, after equilibrium is 
reached, the final wavefunction is given by: 
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For a funnel that uses a single ancilla, the dh\mAt) are the amplitude to start in the state |4i4 100)4*4) 
and transfer an excitation to the bath at mAt. Formally, it can be constructed from the system propagator, 
exp(—iHmAt/Ti), and matrix elements of the system / bath interaction. Using more than one excited ancilla 
results in a two-time dependence for the c. All these functions approach zero for t — > T, due to the irreversible 
loss of amplitude from the funnel states at earlier times. 
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The important point is that the Cn(mAt), for different n and m, uniquely label orthogonal modes of 
the bath. To see the consequences of this, form p from Eq. (Q) by tracing out the bath and environment. 
The populations look like this: 

|V>o)|O0)<M00| |a„| 2 x l4 0) 4°VAi)l4 0) ) + u^<f\mAt)\e^) + • • • | 2 (5) 

I a 

with similar expressions for the other codewords. Because of the earlier conditions on T, the populations 
must be repaired (there is no where else for the populations to go). Thus, the sum in Eq. (||) is one. 
However, from Eq. (Q), it is easy to see that the coherence \ip2)(ipo\ is zero! Using orthogonal ancilla states 
between n=0 and 2 resulting in these codewords becoming entangled with orthogonal bath modes. What 
has happened, is that the pattern of excitation in the bath can be used to determine the probability to be 
in each codeword. Using orthogonal ancilla leaves a separate pattern of excitation behind, which means 
the bath has gained information about the system, and coherence is irreversibly lost jll). Thus, another 
condition for AQEC is that excitation should be symmetrically removed from separate funnels. 
The final condition comes from examining the coherence 

|V>o)|00)(V>i|(00| a* ai x 

£ (4 0) 4 0) (mA0|4 0) ) + uPc£\mAt)\eV) + ■ ■ •) * (u^cf\mAt)\e^) + n«c« (mAt)| e «) + ■ ■ ■) 

(6) 

The sum is the inner product of two vectors, indexed by m, whose elements are environmental wavefunc- 
tions. Each vector individually has a unit norm, so by the Swartz inequality, the sum is one if the inner 
product of each element is maximum. Thus, the final criteria for AQEC is to have J2k c^\mAt)\e^) 
=J2k u q^ c q^ {mAt)\e q k ^) , for each pair of codewords n and q, and at each time mAt. Physically, we are again 
preventing the bath from gaining information about the codewords. In this case, however, the information 
would be transferred by the pattern of environmental entanglements with the bath, instead of the excitation. 

This last requirement is similar to the phase-matching requirement for frequency mixing in non-linear 
optical materials [^Tj. Consider the following contrived example. A qubit suffers an error, a|0) +/3|1) — > 
a 1 2) + /3|3). These four states have frequencies ojq, ui, L02 and U3, respectively. The original coherence 
between |0) and |1) has been transferred to a new set of states. There are no ancilla to this repair; instead, 
relaxation symmetrically drives |2) — > |0) and |3) — > |1) at a steady rate. Therefore, we can write co(t) = 
y/jex-p(—i{u>2t + uJo(T — t)} — ~ft/2) and c\(t) = ^/7exp(— i{uj^t + ui(T — t)} — jt/2). The original coherence 
gains a factor of / Cg(i)ci(i)di = exp(«(wo — oj\)T) x 7/(7 — i(u>o — lo\ — 0J2 + ^3))- When loq — LO2 = u>i — W3, 
the dynamics between the separate funnels is indistinguishable as far as the bath can discern, and a full 
repair results. 

To summarize, the following are sufficient conditions for AQEC, although they may not all be necessary. 
In particular, it is likely that the Markov assumption could be relaxed. (1) The system obeys a Lindblad 
equation of motion, with an evolution superoperator T [111. If cooling occurs on degenerate transitions, they 
must be coupled to orthogonal modes of the bath 24). (2) The eigenstates of H consist of codewords, 
\ip n ), the funnel subspaces associated with each codeword, {|<^n'' ) )}, and the rest. The codewords obey the 



conditions of QECC |L9, 14 , and errors transform codewords only into their associated funnels. If the errors 
are available as joint system / environmental transforms, U, then check whether (ffl \U\4) n ) = (cj) q k) \U\ip q ) 
for all n ^ q, and for some labeling, k, of the funnel states. (3) The codewords form a decoherence free 
subspace with respect to the bath pa]. (4) T does not transfer amplitude between funnels. All funnel 
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populations decay under T into their associated codeword populations. (5) The dynamics under T within 
each codeword-funnel subspace are identical. If ancilla are used, they must be excited symmetrically between 
separate codewords. The last conditions, which are the novel aspect of this approach, are necessary in order 
to repair the coherences of the codewords using dissipation. Alternatively, one could replace the last three 
conditions with criteria on the eigenvalues and left and right eigenstates of T. 



4 Some Numerical Simulations of AQEC. 

4.1 The Single Codeword Model. 

This section provides a better idea of how AQEC works by examining the behavior of a few numerical 
simulations. Recall the example of keeping S in the state |ls). We now see that we could have used (the 
states are ordered as jOsO^), jOsl^), llsO^) and 
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As previously, T has one zero eigenvalue, jlsO^^lsO^j. However, the path by which error correction occurs 



is different: \0s0a ] 



1*1 



\1s0a)- This results in a more efficient repair, as seen by comparison of 



Fig. (H) to Fig. (M). The phase of /z, and the parameters wqi and wio, are irrelevant, but as Ato — ton — luqq 



increases, the first step becomes less efficient. However, it isn't crucial that loqq — 
optimized parameters slow down, but do not halt, the correction process. Note that, 
a spin-flip error at A is removed without ever influencing S. 
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Fig. The linear entropy as a function of time 
during dissipative AQEC for a single stable state. The 
starting state is p = |0s0A)(0s0yi|. The more rapidly 
that H can mix the corrupted state with an excited 
ancilla, and then cool the ancilla, the more rapid the 
repair. Non-zero values of Ao>, or small values of c, 
lead to a slower repair. 



time 



10 



4.2 The Two-Codeword, Spin-hip Correcting Model. 

Let us re-examine the system that protects against spin-flip errors, using three S and two A qubits. The 
codewords are 1 000, 00) and | 111, 00). Parameterize the system H as shown in Fig. (||): 
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Fig. |E|. A parameterized H for a two-codeword 
AQEC with three S and two A. H is block diagonal, 
with the two blocks parameterized as shown above 
(the lines provide a guide for the eye, with a listing of 
the order of the states for the first funnel / codeword 
combination). The blocks must be identical, to within 
a constant offset along the diagonal, so the dynamics 
between the funnels appears indistinguishable. 
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The 7 mix the different error states, the k mix the excited ancilla states, and the fi mix the errors with 
the excited ancilla states. Previous implementations of QECC kept the /i matrix diagonal and the 7, k = 
0, so that separate errors excited orthogonal ancilla states. For AQEC, we must check that each of the six 
eigenstates of the funnels have some non-zero projection along a state with excited ancilla so that population 
is not trapped in a funnel. 

Eq. (^) shows some examples for H. All three sets properly repair spin-flip errors, but set (C), which is 
nearest in spirit to QECC, implements the most rapid repair. For the simulations, only the 14 total codeword 
and funnel states are used in the numerical simulations, so T is 196 x 196 in size. The matrix exponential 
routine of MATLAB |^8| was used to produce exp(Fi). The initial p is found by tracing the environment out 
from the initial error state, |*)|e ) + h,x\^)\ei) + hA^)\ e 2l + £s,x|*)|e 3 ), where |#) = (1/^)1000,00) + 
(exp(z7r/3)/v / 2)|Hl,00). This state allows us to check whether the coherence phase is properly recovered. 
In general, the environmental overlaps (e n \e m ) could be any complex numbers subject to Xl n ( e ™l e ™) = 
and |(e n |e T „)| 2 < (e„|e„)(e 
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Set Eq. (^,A) uses a single ancilla to correct the three independent spin-flip errors. Fig. (|^) shows the 
recovery of the codeword populations and coherences for a spin-flip error at each S, and for a spin-flip error 
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with a set of randomly chosen environmental overlaps, (e n \e m ), as given in Eq. (JsJ) . 
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The H of Eq. (§[A) excites an ancilla only for the first spin-flip error. It repairs the other spin-flip error 
by mixing all the errors together. The numerical simulations shown in Fig. (^) show this process in detail. 
Note that if the all three spin-flips entangle the system with the environment, then the linear entropy of the 
system is initially non-zero. 

0.5i 1 . _ _ . - 0.5 r 





time 

Fig. |^. Repair of the populations (upper left), 
coherences (upper right), and the linear entropy 
(lower), for errors in which each S alone flips, or 
there is a correlated flip of all the spins whose envi- 
ronmental overlaps are given by Eq. (§). All thr ee 
graphs have the same legend. The repair process is 
a linear dynamics, so it also repairs correlated sin- 
gle spin-flips. Populations and coherences of both 
codewords follow identical paths. The parameter 
set is Eq. (||,A), but with the second codeword / 
funnel offset by Aco = 1 in order to show the os- 
cillation of the coherences. The cooling rates were 
ci,c 2 = 1. 

"0 time 10 20 

Thus, AQEC can expell the information about which error occurred at the same time as the error is 
repaired. Eq. (||,B) is another example of this. It mixes together all the errors, and all the excited ancilla 
states. Because it does so symmetrically between the separate codewords, the errors are repaired. Fig. (0) 
shows the populations for all the funnel states after the first spin is flipped. It can be observed that all the 
states are transiently excited: the state vector "swirls around" in each funnel as the repair occurs. 
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Fig. fj]. The populations of the codewords (at 
back), of the three different spin- flip states (mid- 
dle three), and of the three orthogonal excited an- 
cilla states (forward three), during a repair after 
the first spin is flipped. The parameter set is Eq. 
(||,B), with cooling rates Ci,ca = 1. It is permis- 
sible under AQEC to mix together different errors 
during the repair, so long as the dynamics between 
the separate codeword / funnel subspaces is indis- 
tinguishable to the bath. 



state index 



What happens when condition (5) is violated? There are two possibilities: excite the ancilla asymmet- 
rically between the codewords, or have different dynamics between the two codeword / funnel subspaces. 
Fig. (||) shows the first case, for which Eq. (||,A) was used, but modified for the funnel surrounding | 111, 00) 
by setting /in = [in — 1/V%- The cooling rates were c\,C2 — 1, and the error was I\^ x . The partially 
orthogonal ancilla states do not allow Y to return p to a pure state. In fact, if /in were set to zero for 
the second codeword, then there would be no element in V to transfer |000, 10)(111, 01 1 to |000, 00)(111, 00 1 . 
The bath gains information about the system through excitation. The coherence asymptotically approaches 
0.3530 exp(i7r(0. 3333) ), with correct phase but low magnitude. 
O.61 . 
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Fig. A partial repair due to the use of par- 
tially orthogonal excited ancilla states, 1 10) and 
(|10) + |01))/V2, between the two codewords. The 
populations (solid line) are repaired, but the co- 
herence magnitudes are not (dashed line). 




— populations 

— coherences 
linear entropy 



time 5 1 

Fig. |9|. A partial repair resulting from a dis- 
similar dynamics between the two codeword / 
funnel subspaces. The populations are repaired 
(solid line), but the coherence magnitudes are not 
(dashed lines). 



The second possibility is shown in Fig. ([)]). Here, the ancilla are excited symmetrically, but the dynamics 
between the funnels is not equivalent, and coherence is again lost. The parameter set is Eq. (||,C), but /in=2 
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for the second funnel, so the mixing was more rapid. The cooling rates were ci,C2 = 1, and the error was 
I\ tX . Again, the population is repaired, but some coherence is lost. 



5 A Proposed Test System. 



We now give a physically realizable example of an AQEC system that implements Shor's three-qubit, majority 
code against spin-flip errors [jl). ft can not repair phase-flip errors, so it is not suitable for a quantum 
computer, ft serves to illustrate a method by which to find systems suitable for AQEC. It is also encouraging 
that a system can be found without making recourse to exotic interactions. 

Our strategy is to find a system that obeys multiple conservation laws that the errors violate. Consider 
a system with an observable A such that [H, A] = 0, and an error E where [A, E] = E. The simultaneous 
eigenstates of H and A, \e,a), have the property that AE\a) = (a + l)E\a). Thus, choosing codewords 
with a = and 2 gives rise to funnel states with a = 1 and 3. If, in addition, there is a unitary B such 
that B\e,a) = \e,a+ 2), then the funnel states can be mapped onto one another, and their dynamics are 
equivalent. 

Let three spin 1/2 particles be lined up along the z axis, in a zero static magnetic field. They interact by 
point dipolar D nm (I nyX I m ^ x + I n ,ylm,y ~ 21n,zlm,z) and exchange J n m(ln,xlm,x + In,ylm,y + In,z^m,z) terms 
[PI . Dipolar interactions decrease with distance as r -3 , so for equally spaced spins, -Di2=-D23=8.Di3=£, 
where C can be as large as 0.1 cm -1 29 . 

Assuming the dipolar interactions dominate, the level diagram of the spins is given in Fig. @A). The 
spins attempt to mutually align, resulting in ground states of |000) and 1 1 1 1 ) . These are the codewords. A 
spin-flip error, I n ,x, is equivalent to rotating a spin by 7r about the x axis. Since this requires work against 
the dipolar field, dissipation can repair these errors. The funnels come from the conservation of the spin 
angular momentum about the z axis, ^2 n I n ,z, which has eigenvalues denoted as m z . An error changes 
m z — > m z ± 1. The codewords have m z = ±3/2. The funnel surrounding 1 000) (levels A-C of Fig. (|lC|A) ) 
has m z = —1/2, and the funnel surrounding |111) has m z = +1/2 (levels D-F). 
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Fig. 10. (A) The level dia- 



gram for the three spin system. 
The ground states are the code- 
words. States A-C with m z = 
-1/2 form the funnel for |000). 
The dashed lines show the dipole- 
allowed transitions for the first 
funnel, with symmetric transi- 
tions for the second funnel. (B) 
The spectrum of dipole-allowed 
transitions. The starred lines 
represent the funnel to codeword 
transitions that should be cooled. 



Spontaneous emission of a photon with an x polarized B field will symmetrically de-excite the degenerate 
funnels. Thus, we could use photons as the ancilla for this system. There are several advantages to this 
approach. First, it is difficult to selectively cool a single spin, but it is easy to cool an electromagnetic res- 
onator. Second, the validity of the Markov approximation for a damped resonator mode is better understood 
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H [T(|. A y-polarized B field will anti-symmetrically de-excite the funnels, so an error I n ^ y will be "repaired" 
to the phase-flip error I n ^ z . This phase-flip error can not be repaired, because it requires no work to rotate 
a spin about the z axis for this system. 

There are two last difficulties to overcome. The dashed lines of Fig. (|lO|,A) show the dipole-allowed 
transitions for the m z = — 1/2 states. Their strengths are proportional to matrix elements of the operator 
^2, n In,x- Point dipolar interactions alone are not sufficient for AQEC, because (B\ 7 n)X |000) is always 
zero, no matter how the spins are positioned along the z axis. This unwanted symmetry is broken by setting 
J\2 = J\z = and J23 = 0.2£. Actually, any < | J23I < 0.5£ will work. Using the above parameters gives 
rise to the spectrum of Fig. (|Io|,B). The second difficulty is that, besides spontaneous emission of y-polarized 
-B-ficld photons, there are also transitions between the funnels. We wish to cool only the starred transitions 
(at 0.64C, 1.03C, and 2.39C), but not the un-starred ones (at 0.39C, 1.36C and 1.75C). 

This can be achieved by placing the spins at the center of a resonator whose modes are only resonant 
with the starred transitions. Consider a rectangular, conducting cavity of linear dimensions a, b, and d. 
Resonances are indexed as transverse electric (TE mnp , where n + m > 0, p > 0) and transverse magnetic 



(TM„ mp , where n, m > 0, p > 0) modes 1 30 , with frequencies Lu mnp — (m/a) 2 + (n/b) 2 + (p/d) 2 in units 
of cm^ 1 if the cavity lengths are in cm. Each mode produces either a linearly polarized electric or magnetic 
field, or no field, at the center. One can invert the above to find that a resonator with dimensions a=2.32/£, 
6=0.87/C and d=4.28/C, has TE102, TE104, and TE122 modes resonant with the starred transitions. Each 
mode has an x polarized B field at the spins. Another 29 modes exist with u> < 2.5(. Of those that produce 
B fields at the spins, the nearest to a funnel-funnel transition is TE302, which is offset by 0.018£ from the 
C-E transition. A resonator Q ^ 76 is required to suppress emission of this transition. For £ ps 0.1 cm , 
microwave resonators can achieve this goal. This larger £ is also desirable because the cold bath must satisfy 
T <C (hc/k)( » 0.1 K, a not outrageous requirement. 



6 Discussion 

AQEC borrows the same structure for storing information as in QECC, but implements the error correction 
in a different way. In NMR terminology, the qubit of the above system is hidden in the triple quantum 
coherence of the spins. The novel aspect is that dissipation can be used to directly repair not only the 
codeword populations, but also the coherences. The criteria for this is simply summed up by demanding 
that excitation, and environmental entanglements, be expelled from the codewords in a symmetric manner. 

Especially interesting is the possibility that an AQEC qubit exists that can protect against both spin- and 
phase-flip errors. Exchange interactions may prove more useful in this regard. Being isotropic interactions, 
they resist the rotation of a spin about any axis. One difficulty with using only exchange interactions, is 
that dipole-allowed transitions vanish, so a symmetric de-excitation of all the funnels by photons becomes 
problematic. Another open question is how AQEC behaves when it is scaled up to large numbers of code- 
words. It is, however, helpful to contemplate an error correction scheme that requires no additional burden 
to the programmer. 



7 Conclusions 

Conditions are given by which the dissipative evolution of a system, coupled to a cold Markovian bath, can 
be used to implement automatic quantum error correction. The new condition, necessary to repair codeword 
coherences, requires a symmetric de-excitation of separate codewords, and an equivalent dynamics between 
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the different codeword / funnel subspaces. They resemble the conditions of phase-matching in nonlinear 
optics. A test case, that of Shor's majority-code against spin-flip errors [0, is proposed. It utilizes well 
known dipolar and exchange interactions between spins, and dipole-allowed transitions with the modes of a 
resonator. 
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